#!/usr/bin/env python

'''
program plots phase diagram for ferromagnetic triangular transverse
'''

import matplotlib.pyplot
import math
import numpy as np

B = [0,    0.1,   0.2,   0.3,  0.4,   0.5, 0.6,  0.7,  0.8,  0.9,  1,    1.1,  1.2,  1.3,  1.4,  1.5,  1.6,  1.7,  1.8,  1.9,  2,    2.1,   2.2,  2.3,  2.4,  2.5,  2.6,   2.7, 2.8,  2.9,  3,    3.1,  3.2,  3.3,  3.4,  3.5,  3.6,  3.7,  3.8, 3.9, 4, 4.1,  4.2,   4.3,  4.4,   4.5, 4.6,  4.7, 4.76]
T = [3.62, 3.618, 3.612, 3.61, 3.606, 3.6, 3.59, 3.58, 3.57, 3.56, 3.55, 3.53, 3.51, 3.49, 3.47, 3.45, 3.43, 3.41, 3.37, 3.34, 3.31, 3.275, 3.24, 3.20, 3.16, 3.11, 3.055, 3.0, 2.96, 2.91, 2.85, 2.79, 2.72, 2.65, 2.57, 2.48, 2.4,  2.32, 2.2, 2.1, 2, 1.9, 1.77,   1.61, 1.4,  1.2, 0.9,   0.4, 0]


matplotlib.pyplot.text(1, 2, "ferromagnetic",fontsize = 20)
matplotlib.pyplot.text(3, 3.5, "paramagnetic",fontsize = 20)
matplotlib.pyplot.xlabel("B/J")
matplotlib.pyplot.ylabel("T/J")

matplotlib.pyplot.plot(B, T, "*-")

matplotlib.pyplot.show()